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Gamma rays from the annihilation of dark matter particles in the Galactic halo provide a partic- 
ularly promising means of indirectly detecting dark matter. Here, we demonstrate that pronounced 
spectral features at energies near the dark matter particles' mass, which are a generic prediction for 
most models, can significantly improve the sensitivity of gamma-ray telescopes to dark matter sig- 
nals. We derive projected limits on such features (including the traditionally looked-for line signals) 
and show that they can be much more efficient in constraining the nature of dark matter than the 
model-independent broad spectral features expected at lower energies. 

PACS numbers: 95.35-|-d, 95.30.Cq, 95.55.Ka, 29.40.Ka. 



I. INTRODUCTION 

Indirect dark matter (DM) searches aim at seeing an 
excess in cosmic rays from the annihilation or decay 
of DM in the Galactic halo Gamma rays play a 
pronounced role in this respect because they are pro- 
duced rather copiously and directly trace their sources 
as they propagate essentially unperturbed through the 
galaxy. Powerful currently operating telescopes like 
Fermi LAT 0, H.E.S.S. Q, MAGIC [| or VERITAS H 
now start to constrain viable DM models and next gen- 
eration instruments like the planned CTA Q will be able 
to dig quite a bit into the underlying parameter space of 
particle physics models, in a way complementary to both 
direct DM detection and searches at the CERN LHC . 

Very often, indirect searches focus on secondary pho- 
tons from the fragmentation of annihilation products, 
mostly via tt" — > 77. The resulting spectrum is rather 
model-independent and would manifest itself as a broad 
bump-like excess over the expected background at ener- 
gies considerably lower than the DM mass . Convinc- 
ingly claiming a DM detection based on the observation 
of such a feature-less signal will generically be difficult. 

In many models, however, pronounced spectral features 
are expected at the kinematic endpoint = include 
monochromatic gamma-ray lines |8| , shar p st eps or cut- 
offs @, 13 

as well as pronounced bumps The type 
and strength of these features are intricately linked to 
the particle nature of DM; a detection would thus not 
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only allow a convincing discrimination from astrophys- 
ical backgrounds but also to determine important DM 
model parameters (in particular, but not necessarily lim- 
ited to, the value of rn^). So far, only line-signals have 
explicitly been searched for [l^ — despite the fact that 
they are loop-suppressed and thus generically subdomi- 
nant compared to other spectral signatures \U\ . 

Here, we present a general method to search for such 
features and show that these, indeed, help significantly to 
discriminate DM signals from astrophysical backgrounds. 
This allows us to derive very competitive (projected) lim- 
its on both the annihilation rate and nature of DM, which 
we believe will be very useful for DM searches. 



II. METHOD 

The defining aspect of the above-mentioned spectral 
features in the DM-induccd gamma-ray emission is an 
abrupt change of the flux as function of energy; in the 
extreme cases of gamma-ray lines or cut-offs, e.g., the cor- 
responding energy range would simply be given by the en- 
ergy resolution AE/E of the instrument. The basic idea 
that we will adopt here, following traditional gamma-ray 
line searches [12[, is therefore to concentrate the search 
for spectral features on a small sliding energy window 
[Eo,Ei], with Eo<m^<Ei and e = E^/Eq - 0(1-10). 
An important advantage of using small values for e is that 
gamma-ray fluxes with astrophysical origin can often be 
very well described by a simple power-law. In that case, a 
corresponding fit to the data allows an effective determi- 
nation of the background at the statistical limit, greatly 
reducing uncertainties related to astrophysical sources. 

For deriving constraints on spectral features within the 
sliding energy window, we will use a binned profile like- 
lihood method [Hj. To this end, we split [i?o,i?i] in 
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TABLE I: lACT benchmark models that, from top to bottom, 
roughly correspond to the H.E.S.S. 01, the future CTA 0] and 
the proposed DMA 01 telescope characteristics. 




many energy bins Ai?i and define a likelihood function 
i(/x|c) ~ YiiPfj^.{ci) , where {ci) denotes the expected 
(observed) count number in bin i and is the Poisson 
probability distribution with mean /x. Introducing the 
background normalization /3, its spectral slope 7 and the 
normalization of the DM signal a, we have 



10" 



10-' 



10' 



= (dE [dE'VE,E'A,ff{E') 

Sobs J AEi J 



( — - 

dE' 



-(3E' 



1-1 



(1) 



where tohs is the time of observation, Aeff the effective 
area and 'De,E' the energy dispersion of the instrument 
(in the following taken to be Gaussian). Maximizing 
L{iJb\c) for a given data set c results in best-fit values of 
the model parameters a, /3 and 7 within the considered 
window \Eq,Ei\. Upper limits on the signal strength at 
the 95.5% C.L. can then be derived by increasing a from 
its best-fit value until — 21ogL (maximized with respect 
to P and 7) has changed by 4. On the other hand, a 
detection at the 5(t level (neglecting trial factors) could 
be claimed if the best- fit —2 logi values for background- 
only and background-plus-signal fits differ by at least 25. 

Our main assumption here is that the astrophysical 
background locally takes the form of a power law. Obvi- 
ously, this approximation can break down in case of large 
window sizes e, depending on the collected statistics and, 
to first order, on the intrinsic curvature of the background 
flux K = d^ \og{dJBG I dE) / {d\og E)^ : a change in the 
spectral index by A7, e.g., implies roughly |k| ^ A7^/4 
at the transition point; k, could, however, also be affected 
by systematic uncertainties in Aoff. We will derive con- 
straints on the maximally allowed window size £max by 
requiring that these effects do not significantly alter the 
resulting DM limits. 



III. CHOICE OF TARGET AND INSTRUMENT 
SPECIFICATIONS 

For concreteness, we will in the following focus on ob- 
servations of the Galactic center region with Imaging 
Atmospheric Cherenkov Telescopes (lACTs). We con- 
sider the benchmark scenarios summarized in Tab. U 
which roughly correspond to the telescope character- 
istics of the currently operating H.E.S.S. (sj, the fu- 
ture CTA and — as the most optimistic choice for in- 
direct DM searches — the proposed Dark Matter Array 
(DMA) 0. We implement the energy dependence of 
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FIG. 1: Maximal shding energy window SIZG £max as function 
of the window position E. We show for different intrinsic 
background curvatures fCmax ~ 0.1,0.2,0.3 (top to bottom) 
the window sizes above which DM limits are affected by more 
than 50%. The dotted lines show for lACTl £max for which a 
power-law ansatz would still give a good fit to the background. 



the effective area Aes as given in Ref. (Ref. [H) 
for H.E.S.S. (CTA) and take A^f ^ ^ iq . ^cta^ r^j^g 

proton, gamma-ray and electron efficiencies e 



p,7,e 



all 

three scenarios as well as the energy resolution AE/E 
in case of H.E.S.S. and DMA are taken to be energy in- 
dependent; for CTA we adopt results from Ref. [g. We 
will use e-y — e^- = 0.8 throughout and assume that 
the proposed DMA can reject protons with efSciencies 
ep « 10-3. 

For the background, we take into account cosmic-ray 
fluxes of electro ns fisl and protons [l^, [l^ , the diffuse 
gamma-ray flux [ijand the source HESS J 1745-290 ^ 
at (or very close to) the Galactic center. A summary 
and more detailed description can be found in the Ap- 
pendix; there, we also discuss which choice of target re- 
gion Af2 optimizes the signal-to- noise ratio S/Af (see also 
Ref. [lOl). For the Einasto and NEW DM profiles, with 
parameters as in Ref. [2l| , we will adopt a relatively small 
region Afi = 2° x2° around the Galactic center; larger re- 
gions would weaken the signal-to-background ratio, S/B. 
In case of a strong point source-like enhancement of the 
DM signal from the Galactic center (e^ through the 
effect of the super-massive black hole [22|, or adia- 
batic compression [23 - |26j ). it is favorable to focus on 
even much smaller target regions. As an example, we 
consider the case of an adiabatically compressed (AC) 
profile [25I, [2^ for which we choose a target region of 
AO = 0.2° X 0.2°. 

Let us now derive values for the tolerable window size 
£max in presence of maximal background curvatures Kumx- 
Using mock data sets with k = ±Kmax, we compare av- 
erage DM limits (on the DM models introduced below) 
obtained when using the power-law ansatz for the back- 
ground with the limits obtained when incorporating a 
(fixed) curvature k in the background fit. In Fig. [1] we 
display as function of the sliding energy window position 
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TABLE II: DM benchmark models used in our analysis as 
examples for the typical spectral endpoint features to be ex- 
pected in WIMP annihilations. For these particular models, 
we also state the annihilation channel that is most important 
in this context, as well as mass and total annihilation rate for 
thermally produced DM. See text for further details about the 
DM models and Fig. [2] for the corresponding photon spectra. 



E = ^/EqEi, and for different curvatures Kmax, the val- 
ues of Emax above which the DM limits arc affected by 
more than 50%. For comparison, the dotted lines show 
for lACTl the values for e,nax below which the power- 
law fit still appears to be in good agreement with the 
curved background (using as criterion that for at least 
80% of the mock data sets the p-value of the power-law 
fit is larger than 0.05): obviously, a good quality of the 
power-law fit alone does not automatically exclude size- 
able effects on the DM limits. Therefore, a priori as- 
sumptions on K,„ax are indispensable; in our case, we 
employ |k| < Kmax ~ 0.2 - which we checked to be sat- 
isfied for the background we adopt here - to determine 
optimal logarithmic window sizes for lACTl and IACT2 
according to Fig. [1] (for I ACTS, see below). 



IV. DARK MATTER SPECTRAL SIGNATURES 



The DM signal flux from a sky region Ail is given by 
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where (av) is the annihilation rate, dN-^/dE the differen- 
tial number of photons per annihilation, Px{f) the Galac- 
tic DM proflle and s runs over the line-of-sight. For any 
photon spectrum, and a given value of the dark matter 
particles' mass, we can now derive limits on a (aka (crw)) 
by scanning over all possible values of and applying 
the method described in detail in Section |lll As a tech- 
nical remark, we found that the best limits are actually 
obtained by choosing the center of the sliding energy win- 
dow to lie slightly off-set from the kinematic endpoint 
E = of DM spectra at or slightly below which we 
expect to see the features we are looking for. For the 
instrument specifications and background model that we 
adopted here, in particular, the optimal choice turned 
out to be £' = e~^^'^^m^ (not for line signals, however, 
for which we take E = m^). 

In the following, we will discuss three types of typical 
endpoint features that arise from radiative corrections 




X =E/m^ 

FIG. 2: Photon spectra for the DM benchmark models of 
Tab. [TTl Dashed lines show the same spectra, smeared with 
a Gaussian of width Ax/x = 0.1 to give a rough indication 
of how well a detector with such an energy resolution would 
in principle be able to discriminate these models from astro- 
physical (power-law) backgrounds, as well as from each other. 



to the tree-level annihilation process. The most striking 
spectral signature, in terms of a possible discrimination 
from a power-law background, is a gamma-ray line at 
E^f = {E~^ = mj^-[l — m|^^/4m^]), which would result 
from the direct annihilation of DM into 77 (Z7 or H'y) 
[1] . Generically, for thermal cross sections of DM in the 
form of weakly interacting massive particles (WIMPs), 
the annihilation rate is expected to be of the order of 
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cases much stronger line signals are possible [27h29| |. 

As an example for a step-like feature we use the 
gamma-ray spectrum expected from annihilating 
Kaluza-Klein (KK) DM in models of universal extra di- 
mensions [30| . In the minimal version of these models, 
the DM particle is the S'^^, i.e. the first KK excita- 
tion of the weak hypcrchargc gauge boson, and the cor- 
rect rehc density is obtained for m^ci) ~ 1.3 TeV [3l|. 
Its total gamma-ray annihilation spectrum dN/dx (with 
X = E/m^) at high energies is dominated by final state 
radiation (FSR) off lepton final states and turns out to 
be essentially independent of ?77^(i) and other model pa- 
rameters. 

Pronounced bump-like features at i? ~ may arise 
from internal bremsstrahlung (IB) in the annihilation of 
neutralino DM [ll|. While these spectra are in general 
highly model-dependent, we follow here a simplified ap- 
proach by defining two spectral templates dN/dx (which 
we take to be independent of m^) by referring to neu- 
tralino benchmark models introduced in Ref. |ll| . Here, 
BM3 is a typical example for a neutralino in the stau co- 
annihilation region, where photon emission from virtual 
sleptons greatly enhances dN/dx; BM4 refers to a situ- 
ation in which IB from final states dominates. We 
note that the Sommerfeld effect could strongly enhance 
these features, in particular in the case of BM4, in the 
same way as pointed out in Ref. [IJ for line signals. 
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FIG. 3: Thick lines: Expected 2a upper limits on {gv) for selected DM models, DM profiles and observational scenarios; bands 
indicate the variance of these limits. Thin lines: Spectral feature of DM signal has S/B ~ 1% (after convolution with energy 
dispersion). The left panel shows limits on gamma-ray lines, rescaled by a loop-factor of for better comparison. In the 
central panel, the gray band indicates the expected (av) for KK DM, the black part being compatible with the observed relic 
density. In the right panel, we indicate the adopted neutralino benchmark points, and the dotted lines show the projected 5g 
sensitivity. 



In Tab. [Ill wc shortly summarize the properties of the 
DM benchmark models deseribed above, including for 
completeness the actual DM mass and total annihilation 
rate needed to obtain the observed relic density for ther- 
mally produced DM. Note, however, that we essentially 
treat these values as free parameters in our analysis and 
that we are rather interested in the spectral shape of the 
annihilation signal, represented by dN/dx; in Fig. [2] we 
show these spectra for a direct comparison. 



V. LIMITS AND DISCUSSION 

In Fig. [3] we show our results for the expected 2a up- 
per limits (thick lines) on the above DM models as well 
as the variance of these limits among the 300 mock data 
sets that we created for this analysis. We find that in 
particular IB features in the spectrum (right panel) have 
the potential to constrain the annihilation rate at least 
down to values typically expected for thermal production, 
(av) ~ 3 • lO^^^cm^s"^, already for modest assumptions 
about the DM distribution (we verified that NFW and 
Einasto profiles give similar results). This is very com- 
petitive compared to the best current limits from ACTs 
that only rely on secondary photons [32j - though we 
would like to stress that these limits provide rather com- 
plementary information on the DM nature and can thus 
usually not easily be compared. 

For the case of not too strongly pronounced endpoint 
features (like line signals in most models or the step for 
Kaluza-Klein DM), secondary photons will usually be 
more powerful in constraining the total annihilation rate 



{(jv); from the point of view of indirect DM searches, 
however, the detection of the kinematic cutoff will be 
much more interesting than the detection of secondary 
photons since it allows to draw firmer conclusions about 
the DM origin of the signal and even to determine impor- 
tant parameters like m^. For models with very large IB 
contributions like BM3, on the other hand, we find that 
our method provides even stronger limits on (av) than 
what was obtained by the HESS analysis of the Galactic 
center region assuming annihilation into bb (32j . 

In case of an adiabatically compressed profile our lim- 
its could improve by two orders of magnitude, as demon- 
strated for gamma-ray lines in the left panel; under such 
conditions, one could even hope to constrain models with 
very small annihilation rates like BM3 (recall that the an- 
nihilation rate for BM4 is anyway affected by the Som- 
merfeld enhancement (27j and thus likely considerably 
larger than what is shown in Fig. [3]) . As shown in the 
central panel of Fig. [31 the future CTA should be able to 
place limits about one order of magnitude stronger than 
currently possible, and the proposed DMA could further 
improve these by another factor of ten.-'^ 

When probing a specific DM model, the corresponding 



Note that for DMA, as can be seen from Fig. [T] the statistics 
actually become so good that a spectrum with k ~ 0(0.1) curva- 
ture starts to deviate significantly from a power-law background 
already for rather small sliding energy windows. In order to 
obtain reasonable limits, we therefore included k ^ as a free 
parameter in the fit to allow for energy windows somewhat larger 
than shown in Fig. [T] 
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S/B is a good measure for the level on which spectral 
artefacts in the energy reconstruction of the instrument 
must be understood. As can be inferred from Fig. [3] (thin 
lines), most of our derived limits correspond to moderate 
S/B values of at least a few percent (except for IACT3), 
which should be well in reach of current instruments. 

Limits on gamma-ray lines as shown in Fig. [3] are usu- 
ally derived neglecting any secondary gamma-ray com- 
ponent from DM annihilation [l^; this approximation, 
however, breaks down for very small branching ratio into 
lines since part of the secondary component will leak into 
the sliding energy window. Assuming a dominant anni- 
hilation into bb final states, we find that for branching 
ratios into gamma-ray lines smaller than 0(10""*), the 
presence of the secondary flux begins to alter the derived 
gamma-ray line limits significantly. This renders a naive 
application of standard line-search results on DM mod- 
els with generic 0( a!^ ) branching ratios into gamma-ray 
lines questionable |33j . 

The dotted lines in the right panel of Fig. [3] show the 
projected sensitivity to see a 5(T signal in the IACT2 sce- 
nario (neglecting systcmatics and trial factors). Such an 
observation should of course be cross-checked by the non- 
observation of the same signature in control regions with- 
out large DM induced fluxes. A more detailed analysis for 
detectional prospects is beyond the scope of the present 
work and left for a subsequent publication (33j . 



VI. CONCLUSIONS 

Gamma rays from DM annihilation often exhibit pro- 
nounced spectral features near photon energies close to 
the DM particles' mass. Here, we have shown that 
methods from gamma-ray line searches, which greatly 
reduce the uncertainties related to astrophysical back- 
ground fluxes, can successfully be extended to look for 
such spectral features; this provides a probe of the DM 
nature that is complementary to DM searches relying 
only on the rather model-independent spectrum from sec- 
ondary photons. 

While these kind of features may generically be consid- 
ered even more relevant for the detection of DM signals, 
because they would provide rather unambiguous evidence 
for the DM nature of the signal as well as allow to deter- 
mine important parameters like the DM mass, we have 
demonstrated here that including the spectral informa- 
tion may even significantly improve limits on DM signals; 
steps or bump-like IB features can, in fact, be much more 
important in this respect than lines. 

We stress that while we have considered constraints 
for lACT observations of the Galactic center region, the 
presented method is much more general and can be ap- 
plied to both other targets and other instruments; we 
thus expect it to be useful for a wide range of applica- 
tions in indirect DM searches. An obvious extension of 
the approach presented here, finally, is to apply it to the 
detection rather than exclusion of DM signals, as well as 
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Appendix A: Dark matter searches in the Galactic 
Center region 

Imaging Atmospheric Cherenkov Telescopes (lACTs) 
detect gamma rays by measuring the dim Cherenkov 
light produced by electromagnetic showers through the 
atmosphere. Very similar showers are induced by cosmic- 
ray electrons, which hence constitute a practically ir- 
reducible background. Proton-induced hadronic show- 
ers, on the other hand, differ in profile and energy 
density and can currently be rejected with efficien- 
cies tp ~ C'(10~^-10~-^). Due to their large intrin- 
sic fiuxes, charged cosmic rays typically form the ma- 
jor background of lACT observations. For the flux 
of cosmic-ray electrons, we take dJ^- /dEdil ~ 1.17 x 
10"" {E/TeY)~^-^ (GeVcm^ssr)-! above 1 TeV, which 
hardens below 1 TeV to a spectral index of —3.0 [lB| 
(with a transition between the two fluxes that we as- 
sume to be proportional to their generalized mean with 
exponent —2). For the proton flux we take dJp/dEdfl = 

8.73x10-3 {E /TeVy^-'^^ (GeVcm^ ssr)-i [3, which we 
shift to lower energies by a factor of 3 to take into account 
the reduced Gherenkov light output of hadronic showers, 

^rocon. ^ ^true/3 (ggg ^ g p^gf_ 

For observations of the Galactic Center region (GG), 
wc take as further background into account the HESS 
source J1745-290 [H, with dJnESs/dE = 2.3 x 
10-15 (£;/TeV)"^-^^(GeVcm2s)-i. The diffuse photon 
emission measured by H.E.S.S. in a —0.8° < £ < 0.8° and 
\b\ < 0.3° region around the GG is given by dJ^m/dE = 
5.1 X 10-15 (£:/TcV)~^-^^ (GeVcm^s)-! Unknown 
diffuse emission from outside this region will conserva- 
tively be accounted for by upscaling this flux by a factor 
of two within our 2° x 2° target region. We summarize 
all these background contributions in Fig. 21 

The statistical significance of a spectral feature de- 
pends on the signal-to- noise ratio S/N [N ~ ^/B + S) 
inside the target region. The number of expected back- 
ground events B within a target region AJl and energy 
range Ai? is calculated analogously to Eq. ([1]) by replac- 
ing the model flux by the sum of the above background 
fluxes after integrating over Ail. In the same way, the 
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FIG. 4: Summary of adopted background fluxes in a 2° x 2° 
region around the Galactic center, with an effective proton 
flux as it enters the IACT2 scenario. The black solid line 
shows the sum of all background fluxes. 
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FIG. 5: S/M {thick lines) and S/B {thin lines) of DM sig- 
nal inside a circular region around the GC with radius 6, for 
different DM halo proflles. For the Einasto proflle, we also 
show the effect of adiabatic compression (AC). We assume 
S <^ B and an energy threshold of 200 GeV, but for threshold 
energies 100 GeV-5 TeV we flnd similar results. 

number of signal events S follows from the DM annihila- 
tion flux as given in Eq. ([2]). 



In Fig. [5l we show the S/Af (thick lines) of spectral 
features for a circular region around the GC with ra- 
dius 9. We compare results for the standard Einasto 
and NFW DM profiles with parameters as in Ref. [2l| 
{i.e. r^FW ^ 21kpc, rf"- = 20kpc, a = 0.17 and 
= 0.4 GeV cm^^ at Sun's position Rq = 8.5 kpc). As 
can be seen from the figure radii of a few degree are re- 
quired in order to maximize S/Af for these profiles (see 
also Ref. [13). 

For the optimal choice of Ail one should also consider 
the signal-to-background ratio S/B which is shown in 
Fig. [5] for comparison (thin lines). S/B is related to the 
importance of systematic instrumental effects for the sta- 
tistical analysis, i.e. it gives an indication of how well 
artefacts and uncertainties in the reconstructed energy 
spectrum of the instrument must be understood. In most 
of our analysis, we use a relatively small Af2 = 2° x 2° 
region around the GC. As can be seen in Fig.O although 
a larger region could improve S/Af, it also would imply 
a significantly reduced S/B. 

In some cases, e.g. through the effect of the super- 
massive black hole [23] or adiabatic compression |24| - [26| , 
the DM annihilation can be boosted in a region concen- 
trated around the GC, leading to a qualitative change in 
the behavior of S/N' with respect to the above unboosted 
DM profiles. The effect of adiabatic compression is illus- 
trated in Fig. [5] in case of the Einasto profile (Ein.+AC), 
where we excmplarily adopt the adiabatic contraction 
model of Gnedin et al. [2^ together with the best fit 
parameters inferred from the hydrodynamical simulation 
SI of Gustafsson et al. [2^. In this case, the profile inner 
slope steepens to p ^ r~^'^^. For such an enhancement, 
DM self-annihilations start to play a role and constrain 
the halo density to be at most Pmax ~ ^xl [13; 
however, the cutoff radu obtained, r ^ 10"^ kpc, are so 
small that this effect does not influence our results. As 
can be seen from Fig. [5l in the adopted boosted scenario 
it is preferable to consider much smaller target regions 
than in case of the above unboosted profiles; hence, we 
will use a 0.2° x 0.2° region around the GC when calcu- 
lating limits in presence of adiabatic compression. 
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